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Abstract: A systematic study of the different phases of Lorentz-breaking massive gravity 
in a curved background is performed. For tensor and vector modes, the analysis is very close 
to that of Minkowski space. The most interesting results are in the scalar sector where, 
generically, there are two propagating degrees of freedom (DOF). While in maximally 
symmetric spaces ghost-like instabilities are inevitable, they can be avoided in a FRW 
background. The phases with less than two DOF in the scalar sector are also studied. 
Curvature allows an interesting interplay with the mass parameters; in particular, we have 
extended the Higuchi bound of dS to FRW and Lorentz breaking masses. As in dS, when 
the bound is saturated there is no propagating DOF in the scalar sector. In a number of 
phases the smallness of the kinetic terms gives rise to strongly coupled scalar modes at low 
energies. Finally, we have computed the gravitational potentials for point-like sources. In 
the general case we recover the GR predictions at small distances, whereas the modifications 
appear at distances of the order of the characteristic mass scale. In contrast with Minkowski 
space, these corrections may not spoil the linear approximation at large distances. 
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1. Introduction 

Massive gravity has recently received a lot of attention mainly due to its relation to large 
distance modifications of the gravitational force (for recent reviews, see e.g. |[). Even 
if the addition of a Lorentz-invariant mass term to the standard action for the graviton in 
a flat background accomplishes the desired modification, it also implies the appearance of 
new problems, such as the vDVZ discontinuity and the strong coupling of the scalar mode 
of the massive graviton j|, |], ||, It was realized in (see also that some of these 
problems may be softened if the mass term breaks the Lorentz invariance to rotational 
invar iance. 
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It is also known that some of the features of Lorentz breaking massive gravity are 
peculiar to Minkowski space and do not hold in other backgrounds || (see also |L0[ p[|). 
In this note, we will study the behavior of the gravitational perturbations in a curved 
background with a mass term breaking linearized general covariance. 

A caveat to this restriction has to do with the choice of action for the graviton fluc- 
tuations in curved spacetime. Normally, one considers the perturbations of the General 
Relativity (GR) action to second order around a background which solves the equations 
of motion (EOM). This ensures gauge invariance under linearized diffeomorphisms (diff) 
at the quadratic level. In this work, we modify GR by adding a mass term which breaks 
explicitly gauge invariance, but other diff-breaking corrections are possible. In particular, 
also the kinetic term may be modified once one relaxes the constraint of general covari- 
ance. 1 The motivation for considering only mass terms is that we want to focus on large 
distance (infrared) modifications of gravity. Besides, there are known physical examples 
generating this kind of mass terms for gravitational perturbations. 

A first example is a model where the matter sector includes four scalar fields that con- 



dense, breaking spontaneously the symmetry of the background metric |13| (see also [14) 
for related previous work and |15[ for some cosmological implications). In the gauge where 
those scalar modes are frozen (unitary gauge), the spectrum of the perturbations reduces 
to the gravitational modes with a mass term that violates the symmetry of the background. 
In this sense, the scalars are the Goldstones modes of the broken diff invariance. Another 
interesting example is bigravity, where a second rank-2 tensor interacts with the metric 



Qiiv In this case, there are exact flat backgrounds where the metrics do not share the 
whole group of invariance, but preserve a common SO (3). The spectrum of fluctuations 
around these backgrounds includes a Lorentz breaking massive graviton which is a com- 



bination of both metrics fl7| , (see also for some phenomenology and [20, 21] for 
spherically symmetric solutions). 

Inspired by the previous models 2 we will consider the presence in the action of a 
generic mass term (function of the metric) which breaks general covariance. This term 



both allows for FRW backgrounds (see e.g. |24| , 25 1) and generates the LB mass terms 
for the gravitational perturbations. Thus, the analysis of just the gravitational degrees of 
freedom is consistent in this setup where diff-invariance is broken, while in a diff-invariant 
context this is possible only in a de Sitter (dS) space. 3 

The paper is organized as follows: In section ^ we introduce our notations and the 
setup of our investigations. Then we analyze the perturbations with LB terms in curved 
backgrounds for the tensor (section [|), vector (section ||) and scalar modes (section || 
and appendix ^]). In section ^ we study the generalized Newton- like potentials and their 
deviations from GR. We present the conclusions in section P. 



1 Recently, there has been some interest in the modifying the kinetic structure of GR as a way to improve 
its UV behavior Jl|]. 

2 Other related models include theories with extra dimensions 0] and theories with condensing vector 
fields ||. 

3 We leave the analysis of gravitational perturbations coupled to additional fields in a FRW background 
for a forthcoming publication 
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2. Action, background and perturbations 



Our starting point is the Einstein-Hilbert (EH) Lagrangian with the addition of mass 
terms for the gravitational perturbations, breaking general covariance. In the flat limit 
these terms also break Lorentz invariance, and we will refer to them as Lorentz-breaking 
(LB) terms. 

This setup describes, at quadratic level, infrared modifications of gravity where only 
gravitational degrees of freedom are present. At full nonlinear level these deformations of 
the EH theory may be parametrized by adding to the lagrangian a nonderivative function 
of the metric components, breaking general covariance: 



S = Jd 4 x^M 2 P [R-2F( gilv )\ . (2.1) 



General covariance can be restored [||, |8| by introducing extra (Stiickelberg) fields, 4 which 
in the equivalent of the unitary gauge yield the form (|2.l| ) . 

It is clear that the term F will contribute to the background EOM, and exact solutions 
are known for certain F functions. For example, when F = 2A = constant, the homoge- 
neous and isotropic (FRW) background will be maximally symmetric and the theory will be 
gauge (diff) invariant. For F ^ constant, FRW solutions can be found, which can modify 
the standard cosmological solutions of GR. For certain classes of F, solutions that exhibit 
late-time cosmic acceleration were studied in [jl3| (see also pq| ). 

Accordingly, we assume that the dynamics of modified gravity admits a spatially flat 
isotropic and homogeneous background (FRW henceforth) 

9fj,v = a(r]) 2 r]^ u with rj^ = diag(-l, 1, 1, 1) , (2.2) 

where r\ is the conformal time. We will use TC(r]) = a'/a and Hirj) = a'/a 2 , where ' is the 
derivative with respect to 7] (therefore aH' = TL' — H. 2 ). 
We define the metric perturbations as 

9 tiv = a 2 (rjfu, + h^y) . (2.3) 
The second-order expansion of S can then be written as: 

S = 4« + 4 2 1 > (2-4) 

where 5 

S GR = ld 4 x^Mj,(R-6H 2 ) (2.5) 



and the term F gives rise also to LB masses for h. Assuming rotations are preserved, these 
can be parametrized as 



,(2) _ M 2 



ml h 2 m +2m\ ^-(rrft-AH' 'a" 1 )/if j +(m^-2i7 / 'a~ 1 )/i? i -2m| h 00 h 

(2.6) 



4 This implies the addition of (at most) four scalar fields, but invariant actions can also be found by 
adding vector or tensor fields. 

5 We stress that in this expression, H depends only on the background. 



-3- 



Here spatial indices are contracted with 5ij, and rrii = m^rf) represent effective time 
dependent masses. The terms proportional to H' in ([0|) are conveniently chosen to cancel 



similar contributions coming from the expansion of (2.5) in backgrounds different from dS. 
Notice that the parametrization in ( |2.6| ) is completely general as the mass parameters are 
arbitrary functions of the conformal time. 

Diff gauge invariance is restored taking the limit m, = 0, H' = and it corresponds to 
the case F = A =constant. On the other hand, for FRW background and vanishing masses 
the action is invariant only under longitudinal spatial diffs. This is a consequence of the 
fact that a generic FRW background is never a consistent background for GR without 
matter. A non maximally symmetric background breaks the time diffs, and accordingly in 
the limit of vanishing masses one recovers the invariance under spatial diffs. 

In the Lorentz-invariant case the masses can be expressed in terms of two parameters 

a, 

mo = a + (3 , m\ = —a , m\ — \aT x H' = — a , m| = /3, m 3 — 2cT x H' = j3 , (2.7) 

and the mass term ( |2.6| ) can be written in terms of contractions of h„ v with g^ v '. The 
Fierz-Pauli (FP) choice, free of ghosts in flat space, corresponds to a + (3 = 0. In curved 
space also the "non-Fierz-Pauli" case a + / can be free of ghosts (see section |6|). 

The setup introduced here is suitable to describe a rather general class of massive 
gravity theories, and exhibits a rich set of phases, depending on the masses m, and H{rf). 
In general we can have the following scenarios: 



• The F term in the action (2.1) does not affect neither the background (dS limit) nor 
the propagation of the perturbations. This is the case for H' = and m; = and is 
realized, e.g., when F = A. 

• Only the perturbations are modified. This corresponds to H' = and rrii ^ 0. 
Writing F(g fiu ) = A + /(g^v), this happens when the scale of A is much larger than 
the scale related to /. 

• Only the background is modified. This happens for H' ^ and rrii = 0. As we will 
see, it can also be realized in less trivial situations. 

• Finally, in general both the background and the perturbations are modified by F. 

In order to study the dynamics of the perturbations, it is convenient to decompose the 
metric fluctuations as irreducible representations of the rotation group: 6 

^oo = ip , 

hoi = Ui + diV , diUi = 0, (2.8) 

hj = Xij + diSj + djSi + didjcr + r, <9;S; = djXij = ^ijXij = • 



Here, we follow the notation of 
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From those fields one can define two scalar and one vector gauge invariant quantities 
* = r + H (2v - a') , $ = ib - 2v' + a" - H (2v - a') , 

(2.9) 

W i = u i -s i , 

while the transverse-traceless spin two field Xij 1S already gauge invariant. It is also con- 
venient to define the field £ = a/ A. 

We will couple the gravitational fields to a conserved 7 energy- momentum tensor T^, 

S T = - J d 4 z a 2 V, rT^T^ = ~ f a 2 ( Xij Ty + $T 00 - 2 T 0i W % + * T«) , (2.10) 

with g au V ' a T^ v = 0, where V is the covariant derivative associated to the background 
metric. The field is the generalization of Newtonian potential around the source in the 
linearized approximation. For FRW, the EMT conservation is equivalent to 

U = diT 0i - H (Too + TijSij) , djTij =T[ + 2H T l0 = a' 2 (a 2 T l0 y . (2.11) 
3. Stable perturbations in curved backgrounds? 

Our goal is to study the dynamics of perturbations in curved backgrounds and determine 
when one can get a theory free of instabilities. These instabilities can be of ghost or Jeans 
type. The ghost-like instabilities are related to an infinite phase-space volume. If they 
are present, the decay rate of the perturbative vacuum will be infinite unless a cutoff is 



introduced in the theory []28|| (see also [29]). In the Lorentz breaking case, the different 
masses provide a natural energy scale to place the cutoff. If we admit a hierarchy inside 
the mass scales we may freeze the non-stable degrees of freedom, while still keeping some 
of the masses below the cutoff. In an expanding universe, there is also another important 
dimensional parameter, H. We will focus on modifications such that at least some of the 
mass scales are inside the horizon scale mi S> H. In this case, there is a natural hierarchy 
inside the set of dimensional parameters, which allows to define a large momentum cutoff 
A c keeping the masses small, 

|A|<W=y~A c °, (3.1) 



with a > 0. Even if the addition of a cutoff may unveil phenomenologically acceptable 
phases, to keep the discussion simple we will consider theories free from ghost-like instabil- 
ities in the quadratic Lagrangian (see however [||, |i"7|), 

Jeans-like instabilities can be also present. By this we mean instabilities that appear 
in a certain finite range of momenta. They are the signature of the growth with time of 
the perturbation at certain scales, and may even be interesting phenomenologically as a 
contribution to the clustering of matter at large distances (see, e.g. (3(J). Furthermore, in 
an expanding universe, they may be settled beyond the horizon, where they are presumably 
frozen. Again, we leave the study of this possibility for future research p6| , and concentrate 
on lagrangians with a stable spectrum. 



7 EMT conservation is not strictly required in massive gravity. A study of non-conserved EMT in FP 
can be found in 12711 • 
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It is also important to recall that in a FRW universe, energy is not a conserved quantity, 
and its positivity does not guarantee stability. Nevertheless, for scales smaller than the 
horizon, we can still use the positivity of the energy associated to the conformal time as 



a necessary requirement for stability (see also [31, |3^[). As we will focus on these length 
scales, we will not discuss any global issue. 

Finally, there are two concerns for massive gravity beyond the linear theory The 
first one is the strong coupling that emerges when one or more propagating states have 
their kinetic terms suppressed by a small parameter. In this case, the range of validity of 



the linear theory is drastically reduced. Furthermore, if we consider the action (2.1) as 
an effective action with a cutoff, one expects the contributions of higher order operators 
to become important at much lower energies than the initial cutoff scale. To study this 



behaviour one should analyze the scaling of relevant interaction terms [33] which is beyond 
the scope of this paper. Here we just point out that when Lorentz invariance is violated, 
the strong coupling cutoff can be present in energy and/or momentum independently. We 
will accordingly speak of time and space cutoff A t , A s , by making canonical the relative 
quadratic terms in the action. 

Besides, when the functions mi satisfy certain conditions, there is a reduction of the 
phase-space, i.e. not all the six DOF of the gravitational perturbations propagate. It turns 
out that in Minkowski those are the only ghost-free possibilities ||. In general, unless there 
exists a symmetry that enforces them (see for example the case of bigravity [18|, or the case 



described in jl3|]) these conditions are only satisfied for the quadratic Lagrangian in very 
finely tuned backgrounds. This means that the analysis is very sensitive to small changes 
in the background and probably to the interaction terms and higher order operators Jl|, 
|, |l|. In the following we show how the generalization of LB massive gravity to curved 
backgrounds is useful to circumvent both concerns. 

In the next sections we analyze the spectrum of tensor, vector and scalar perturbations, 
that at linearized level are not coupled, by SO(3) symmetry. 

4. Tensor Modes 



The action for the tensor perturbations is 



il^d^XijduXij ~ a 2 mlxijXij , (4.1 



from which the EOM read 



V^d^Xij ~ Wx'ij ~ o?m\xij = . (4.2) 
The absence of tachyonic instabilities requires 

ml>0. (4.3) 
One can also readily check from ( |4.1[) that there are no ghost or gradient instabilities. 8 



In an expanding universe, the friction term appearing in (4.5) implies that the perturbation is frozen 
at large distance. When imposing mi > 0, we are assuming that the mass scale is well inside the horizon 
and unless otherwise stated, we will assume this to be the case. 
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5. Vector Modes 



Extracting the vector part from (2.4) we get 

M 2 f ( i 

S (v) = _p / d 4 x 2| _ ( n . _ gQAfa - s'i) + a 2 [ml um + m 2 2Sj A Sj ] } . (5.1) 

The field m is not dynamical and it can be integrated out through its equation of motion, 

A(m - s^) - a 2 m\ui = 0, (5.2) 

to yield 

7 



M ? I d*x a 4 



A 



js'i + m 2 SiAs t 



(5.3) 



A — a 2 m 2 ' 

Therefore, the dispersion relation of the vector field Si breaks Lorentz invariance at any 
scale (provided mi ^ m-i). The action is free from instabilities for 



m 1 > and 



m\ > 0, 



(5.4) 



in complete analogy to what happens in Minkowski space. 

The case m\ = is particularly interesting, as it implies the cancellation of the time- 
derivative term in (|5.3j), so that there is no propagating vector mode. As we will see in 
section |6.3| , this case is also important for the scalar sector. 

The canonically normalized field s\ can be defined by a rescaling: 



A v (r)) 2 = a 2 m\M P A 



A 



A — a 2 m 2 



with the action 



= J / d 4 x { WW + ^sfAsI - a 2 m 2 + A 2 V (A~ 2 ) 



m\ 



-2\ll 



(5.5) 



(5.6) 



Therefore the canonical field has a LB dispersion relation and a time dependent mass. 

From ( |5.5| ) we can also read the naive temporal strong coupling scale of the vectors, 
that is momentum (and time) dependent. At large momenta, |A| > m 2 a 2 , we expect 
the physical strong coupling scale to be given by A t = A v /a ~ ^Jm\Mp, that is also the 
expected cutoff for a gauge theory explicitly broken by a mass term m. This implies that 
the theory can be trusted only if the horizon scale does not exceed the cutoff, H < \fm\Mp. 
In section |8| we will comment on some physical consequences of this bound. A similar spatial 
strong coupling scale can be defined by making canonical the spatial gradients. 



6. Scalar Modes 

The scalar sector of the theory is the most interesting one. In flat space with the FP 
Lorentz invariant mass term it is a scalar mode which has the lowest strong coupling scale, 
and it is this sector that shows crucial differences when Lorentz symmetry is violated or 
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spacetime is curved [[?], [|. As we will see later, it is also here that the difference between 
the maximally symmetric spacetimes and generic FRW spaces arises. The analysis will 
show that generically there are two scalar degrees of freedom, while the only possibilities 
with less degrees of freedom are m\ = or mo = 0. 

The scalar part of ( |2.4| ) can be written as (modulo total derivatives) 

s (s) = yd 4 xa 2 {-6(T' + Hip) 2 + 2(2^ - r)Ar + 4(r' + H^)A{2v - a') 

m\^ 2 - 2m\vAv - m\{pA 2 a + 2rAcr + 3r 2 ) (6.1) 



+a 2 



+ml(Aa + 3r) 2 - 2raf^(A<7 + 3r) } , 



where the H , if' terms have canceled as promised. 

In the de Sitter background (dS), when all masses are set to zero, the action reduces 
to the first line of (^O]) and it is gauge invariant. 9 As previously remarked, for a FRW 
background and vanishing masses the action is invariant only under longitudinal spatial 
diffs (only a is undetermined). 

From ( |6.1| ) it is clear that ip and v are Lagrange multipliers enforcing the following 
constraints 

m\m\ (S + 3r) a 3 + 2Hm\ (£' + 3r') a 2 - 2Am\ra - 8H At' 



8H 2 aA + (m 2 . - 6H 2 ) m?a 3 

(6.2) 



2Ha 2 m\ (S + 3r) + 4F 2 aS' + 2m 2 ar' - AH At 



8H 2 aA + (mg - 6F 2 ) mfa 3 

Notice that the behavior of these fields in FRW is qualitatively different from Minkowski 
space. In particular, whereas in flat space, the cases mo = and mi = are singular 
and must be treated separately, in curved spacetime, ip and v are always determined by 
equations Q6.2j ). After integrating out v and ip we are left with a Lagrangian for ip = (£, r)*: 

£ E)T = i//C^+ ^BifZ-^Af, (6.3) 



where 



K _ -M 2 a 2 /2H 2 a 2 m 2 a 2 m 2 m 2 \ 

8H 2 A + (m 2 - 6H 2 ) m 2 a 2 \ a 2 m 2 m 2 m 2 (3a 2 m 2 - 4A) J ' 



Mj, AH (mj - 2m 




8if 2 A + (mg - 6i? 2 ) m 2 a 

We will first study the dynamics of ( |6.3D through the Hamiltonian, for which the explicit 
expression of the matrix A is not needed 10 . The conjugate momenta n are 

dC 



9 The first line of ( |6.l[ ) differs from the standard action of the graviton in a FRW background by a term 
proportional to (Tif — TL 2 )ip 2 which cancels in dS (cf. Jji|). 

10 As the form of this matrix in the general case is quite cumbersome and not particularly illuminating, 
we will not write it explicitly in this work. 
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Thus, two DOF will propagate when the matrix K, is non-degenerate, i.e. when 



det | |/C| | oc mo m\ / . 



(6.7) 



In this case one can express the velocities in terms of momenta, and the resulting 
Hamiltonian is: 



H s , r = ^JC-W + ^Mtp, M = (A + BK- l B) , 



(6.8) 



with a rather simple kinetic term: 



KT 



M 2 a 2 



\ 



4A 

a 2 m? 



2iP 

m 2 Q J 



(6.9) 



The theory is free of ghosts when the kinetic energy matrix K, 1 is positive definite, that 
translates into the following conditions: 



m\ > . 



< ml < 6H 2 . 



(6.10) 



Therefore contrary to the flat space case, we can still have a well defined kinetic term with 
two propagating degrees of freedom. In fact a window for ml opens up, and this allows 
even for a "non-FP" Lorentz-invariant mass term free of ghosts (and vDVZ discontinuity, 
see section [7]). 11 

It is also instructive to look at the no-ghost conditions in the low and high momentum 
regimes. We find, 



no ghost 



at large momenta: m\ > , ml > 

at small momenta: m\ > , < ml < 6H 2 . 



(6.11) 



Therefore a nonzero curvature allows the scenario where the theory is free of ghosts in the 
ultraviolet but there is one ghost mode at large wavelengths; this happens for m\ > and 
ml > 6H 2 . Such a ghost mode at very large distances would not necessarily render the 
theory phenomenologically sick, but would indicate a large scale instability of backgrounds 
with curvature smaller than ml/ 6, including the limiting case of Minkowski (Jeans-like 
instability, in the language of the section ||) . 

From ( |6.9D one can find when the scalar sector suffers from strong coupling due to a 
small kinetic term. As happens for the vector modes, one of the strong coupling scales is 
related to the smallness of mi, whereas the other one depends on the ratio ml/H 2 . When 
this ratio is not small (and compatible with the ghost-free condition ( |6.10| )), both the scalar 
and vector sector become strongly coupled at the same time scale ~ \fm\Mp. 



11 Recently, non-Fierz-Pauli lagrangians with scale dependent masses were also considered in pq ]. Notice, 
though, that in that case Lorentz invariance made the masses depend on both space and time, whereas in 
this work we are dealing only with time dependent masses. 
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The analysis of the positivity of the "mass" term Ai is rather cumbersome and we will 
consider just the high momentum limit (larger than the rest of the scales: rrii, H and H'). 
In this case, requiring that the mass matrix in ( |6.8| ) is positive definite gives 



ml 



m\ < 



to; 



2ml) 2 



16to-q 



H'a 



< 



TOT 



TOT 



2\ 2 



2ml) 



16to 2 



(6.12) 



(where we have used m\ > 0, TOq > and H 2 > 0.) When the previous conditions are 
satisfied there is no gradient instability at small distances. Notice that the r.h.s. in the last 
condition is always negative, meaning that only a FRW background with an expanding 
horizon can be stable. Besides, one can easily check that the previous conditions are 
inconsistent in the Lorentz-invariant case. 

To summarize, in the non degenerate case of too,i ^ 0, we found two DOF where 

• there is no ghost provided m\ > 0, 6H 2 > m\ > 

• there is no gradient instability when ( |6.12 ) are satisfied (so H' is negative). 



The difference that we found with the maximally symmetric case, where there is necessarily 
a gradient instability, implies the presence of a spatial strong coupling problem in this 
limit. In fact in approaching the dS background the spectrum of the Hamiltonian must 
pass through the case in which one of the modes is frozen, because the determinant of Ai 
vanishes and accordingly the "spatial" part of its dispersion relation will vanish. 12 

In the degenerate cases, when mo or m\ vanish, there are less DOF and a separate 
analysis is given in the following sections. The too = case is related to the Fierz-Pauli 
case p6| , whereas the case m\ = appears naturally in the ghost condensate and bigravity 
theories gj| 

6.1 The phase too = 

For too = 0, the field r is an auxiliary field as one can check from the action (|6.3| ). Even if 
there is only one remaining DOF, the general treatment is quite involved and it is presented 
in appendix |A]. In this section we will just state the results and study some particular cases. 



The EOM for r yield the constraint ( A.l ), which once once substituted in the action 
gives a (quite complicated) effective Lagrangian for X. Its kinetic part is 

K, 3a 3 ml [a{m\) 2 + 2(am 2 H 2 - H(m\)' + m\H')] - 4\am\(m\ - m\) + mjH']A 



a 4 M 2 oto 2 (2A - 3a 2 TOl) 2 - 2(4A - 3a 2 m 2 )[3a 2 H(aHm 2 l - (to 2 )') - (2A - 3a 2 ml) H'} ' 

(6-13) 

where to 2 = 3(to 2 j — to 2 ) — to 2 ,. The positivity of the kinetic energy (no ghost) for large 
momenta gives 

am\(m\ — m-1) + m 2 H' 



amj + 



> 0. 



3.14) 



12 See || for a discussion of the modifications to these dispersion relations coming from higher order 
operators. 



- 10 - 



On the other hand at small momenta the kinetic term reduces to 

£|a=o = a 2 Mj,/3, 

which is, remarkably, always positive. 

One can show (cf. appendix [A]) that K, is positive also at any intermediate momenta 



in the variable A provided that, in addition to ( 6.14 ), one has 



/ a(m\ + 2H 2 m 2 ,) + 2m 2 d H' - 2H[mT)' \ 

to 2 >0, 27"^ 2T 2rr/ >0- (6-15) 

\ am\(m\ — m\) + mfH' I 

When these conditions are saturated we are led to a case with a vanishing kinetic term, as 
discussed below (see ( 6.24D ). Due to its analogy with the special case discussed in [32], we 



will refer to this case as partially massless. 

The condition ( |6.15 ) refers to modes at large distances (eventually outside the horizon) 



and is not present in the Minkowski spacetime. In dS, taking the Lorentz-invariant FP 
limit (0) with m 2 = (3 = —a, the previous conditions reduce to the Higuchi bound f37|. 



2H 2 < to 2 . (6.16) 

Contrary to this case, the LB mass terms allow for a unitary massless limit. In fact, if the 
mass is an appropriate function of the conformal time, this limit can be free from ghosts 



also in the Lorentz-invariant case (see section 3.2.3). 

From the same kinetic term we can also estimate the strong coupling scale of the field 
E, since the canonical field S c is defined at high momentum by the rescaling 



A V A MpTOl 



am? 



v^A 
where x = m\jm\. 

Concerning the potential term M, it can be written as 



(x(l-x) + H' 
ami + AH' 



(6.17) 



m 



b 2 + cA + dA 2 + eA 3 + (to| - m|)A 4 



M = — ■ -= ^ (6.18) 

T 

where 6, c, d, e are functions of TOj and Ti., whereas q is a second order polynomial in A. 

The absence of gradient instabilities, equivalent to the positivity of M, requires in the 
ultraviolet and infrared regimes the following simple conditions: 

{at large momenta to 2 . > m\ 
3 (6-19) 
at small momenta to 2 > . 

We see that for to 2 . > m\ the potential is free from gradient instabilities that would be 
as dangerous as ghost instabilities as they would imply an infinitely fast instability p9[ . 
Notice also that at zero momentum, the condition to. 2 > required for the stability of 
tensors and vectors, enforces positivity of the potential. This implies, together with the 
stability of the kinetic term, that at small momentum the theory is always stable. 

At intermediate scales, the analysis becomes very technical, and a method to check for 
the positive definiteness is presented in the appendix |A[ 
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6.2 Particular cases with mo = 

Some interesting subcases of the tjiq = dynamics can be found looking at the numerator 
and denominator of eq ( 6,13; ). When the denominator vanishes, the field r disappear from 



the EOM and the constraint ( |A,1| ) does not hold anymore. The analysis of this situation 
is presented in section 6.2.l[ 



When the mass parameters are fine tuned in such a way that the numerator of ( 6.15Q 



vanishes, the £ field does not propagate and it becomes an auxiliary field. This possibility 
is examined in section |6.2.2 , where we also show its relation to a gauge invariance related 



to conformal invariance. Finally, the Fierz-Pauli Lagrangian in dS is a subcase of the phase 
niQ = where this fine tuning can occur. We study this possibility in section 6.2,3| . 



6.2.1 Time diffeomorphisms in FRW 



From the constraint ( A.l ) for r as a function of S, we see that it is singular for specific 



values of the masses. This happens when 

T2V 



m 2 = 2m\ = -Aa^H' , ml = a ^ H \ TT 2H . (6.20) 

a 2 H 



(If only the first condition holds, the constraint for r ( |A.l| ) reduces to r = — S/3.) 

In this case, and away from dS (we are assuming m\ ^ 0), the field r does not appear 
at all in the action. This corresponds to a restoration of the gauge symmetry corresponding 
to time diffeomorphisms, 5{2v — a') = 2£o- 

The final Lagrangian in terms of £ is then given by 



Mia 2 ( a H' 1[ ,. . 2>{aH') 2 + (2aH' + H"/H)A 



A + 3aH> 3 



mi + A- 



(A + 3aH 



l\2 



} . (6.21) 



Notice that in Minkowski this phase has mj = m\ = m^ = and features an enhanced 
gauge invariance mentioned in [18]. In contrast to the case of Minkowski, and to the case 



of dS, in FRW background the field £ propagates. It is also clear from ( |6.21 ) that the 



kinetic term is positive definite provided that H' < 0. Concerning the potential term, 
let us consider scales well inside the horizon. The requirement of positive energy at these 
scales gives the condition 

>-(2aH' + H"/H), (6.22) 



in 



When inequality (|6.22| ) is exactly saturated, the scalar degree of freedom has vanishing 
speed in the at high momenta. Its dispersion relation is then to 2 ~ const. + 1/A. 

Finally, let us note that if for H' < the scalar mode is well behaved, the limit of 
vanishing H' leads to a vanishing kinetic term and thus to strong coupling once interactions 
are taken into account. The resulting time strong coupling scale can be estimated as 



A s = a MpJ - aH ' — . (6.23) 
V A + 3aH> 1 ; 

that is clearly more dangerous at short distances |A| S> \aH'\, where it may become sensibly 
lower than Mp. 
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6.2.2 Partially massless 

Another particular case appears when, after integrating out r, the kinetic term of £ cancels. 
This happens when the inequalities (|6.14j) and ( |6.15|) are saturated, 

4 

m\ = — a I riA — , a(m\ + 2H 2 m 2 u ) + 2m\H' - 2H{m\)' = 0. (6.24) 
am\ + H ' p 

In the Fierz-Pauli limit this expression reduces to the partially massless case of de Sitter 
space(cf. |32|) and corresponds to a situation without propagating scalar degrees of freedom. 
In the Lorentz invariant case with constant masses in de Sitter space, this fact is related 
to a conformal invariance [38]. In the most general case, one can prove that the system is 
invariant under the transformation 

5if> = -2(? + HO + <t>t, 5v = -£ + (', 8r = ZH€ + </>., 5a = 2(, (6.25) 

with 

t _ al(aml( + H'( + H(') _ , 



H(aml + 2H') ' Ys ~ ^ 



"4 

.2 



Cm 2 (am 4 (m2 - AH 2 ) + 2m 4 H' - AHm' A ) 



(6.26) 



2aH 2 

only when the extra condition 

a[m\ + 2B 2 (2m\ - 3ml)] + 2\m\H' - H(ml)'\ = (6.27) 
is satisfied. The previous condition implies the cancellation of the potential part once 



( 6.24 ) is satisfied. In the Lorentz-invariant limit with constant masses and dS background 



(6.27) is always satisfied when (6.24) holds. Notice also that the existence of this sort of 
scale invariance is general even if the kinetic term is not invariant under diff away from de 
Sitter. 

6.2.3 Lorentz-invariant FP limit with time dependent masses 

In the Fierz-Pauli limit ( Q2.7D with m 2 = (3 = —a) the mode S propagates. However, the 
conditions ( |6.24 ) can be still be satisfied in dS (and only for this background) provided 



that m satisfies the differential equation 

mm' = {a 2 m 2 - 2H 2 )m . (6.28) 
This equation can be integrated to yield 



m2w = ^fT(2^^fj^)' (6 - 29) 

where mi is the value of the mass at the time corresponding to a{rf) = 1. The resulting 
mass runs from m 2 to 2H 2 when a runs from to 1. Notice that choosing the initial 
conditions corresponding to a constant mass, m 2 = 2H 2 , we recover the partially massless 
case discussed in |J2| . A similar situation could be studied for the non-Fierz-Pauli (Lorentz- 
invariant) case (mi = m 2 / m^ = m.4, mo = m^ — m<i). 
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6.3 The phase mi = 

The case m\ = is particularly interesting in the Minkowski background, as only the 
tensor modes propagate. As we will show, there is a corresponding effect in dS, while one 
scalar mode starts to propagate in a FRW background. When m\ = the fields £ is not 
dynamical as one can check in action (|6.3|). Accordingly, its EOM is 



W(m!-m§)E 



m 2 r' 



r H{m\ — 3m|)r . 



(6.30) 



Notice that again the Minkowski space limit TC = is peculiar and the degree of freedom 
associated to r is not present 13 . In curved space, generically £ is determined by ( |6,30 ) and 
when it is substituted back in the action, after integration by parts, yields the Lagrangian 



C 



M 2 a 2 



mz 



>2 



H 2 I 2(m 



777-5 



H' 



A 



m 2 2 [H 2 



(ml 



3771? + 3m 2 ] 



m 2 H'a] 



m% 



H[m\{ 



mZ[m 



2y 



771 a 771 



2\A 



771$ 



+ 777,2 



(6.31) 



777-5 



777-5 



777 



2Vl 



( 



777-5 



7775 



where m^ = m\(m\ — 777-3) + 7774. Prom the previous expression we discover that in dS, 
the phase m\ = has no propagating degrees of freedom (in the sense that the action 
is A independent so that there is no dynamics in x space), even if, in comparison to 
the Minkowski case, the scalar sector has a kinetic term from which we expect a ghost 
condensate like dispersion relation coming from higher derivatives ||] . Besides, the potential 
strong coupling scales A s and At are easily read out from the previous expression. 

Thus, in general the phase m\ = is quite rich, and particularly simple. Ghostlike 
instabilities are avoided imposing 777^ (ml — 777-3) — 0- To get r id of gradient instabilities in 
this case, it is enough to impose H' < 0, whereas the tachyon free condition can also be read 
from ( 6.31| ). For the case with constant masses, it reduces to ml\H 2 (m\ — 2>m\ + 3t77|) — 
mjH'a] > 0. 



6.4 Particular cases with m\ = 

A direct inspection of ( 6.31 ) and ( |6.30 ) shows some interesting subcases for the mass 
parameters. First, when the r.h.s. of eq. ( 6.3C| ) cancels, this equation is no longer a 



constraint for S. Besides, for 777^ = 0, the kinetic term for r cancels in the action. We 
devote the rest of this section to the analysis of these possibilities. 



6.4.1 The case mo 



7775 



When ml = m 3 the kinetic term of r is zero. In this case r is non-dynamical and can be 
eliminated from the action. The only degree of freedom now is £ with a Lagrangian 

a 4 Mp (677727774 — 97774 + 4777q 777,2 



777^ 



(2m 2 



3m 2 ) 



(2m^ml 



3m|)' 



777 



v2 



2 (27772 



37772) H 2 



2 v^2 
?77 9 2j 



(6.32) 



Again, this mode has no dynamics in space. From direct inspection we can derive the 
strong coupling scale, and the region of parameters where this mode disappears. 



3 Also, the case = m-j, should be treated differently. 
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6.4.2 The case m v = 

Finally, for m, = we are back to a situation without scalar propagating degrees of freedom 
but still with a potential part at the linear level. In Minkowski also this part vanishes and 
the field r is not determined (indeed, there is an additional gauge invariance). In dS, this 
happens when 

_ n - 3m| + 3m 2 ) ~ ( m lY _ m 2 3 (m 2 )' - m 2 (m 2 )' 

m V ~ U ' 2 — 27 2 2\ ' 1^0. OOJ 

and outside this region of the parameter space, the EOM gives r = 0. 



7. Coupling to matter and vDVZ discontinuity 

Though the vDVZ discontinuity is one of the main phenomenological difficulties of FP 
massive gravity in flat space, it is known that it may be circumvented in curved backgrounds 
or when one considers Lorentz violating mass terms [j?]]. For AdS or dS, the vDVZ 
discontinuity is avoided by hiding the effects of the mass at distances larger than the 
horizon, and as a consequence there is no modification of gravity at scales smaller than 
the Hubble radius. In this section we will see that some of the massive gravity phases we 
have studied allow for a modification of gravity at scales shorter than the horizon scale 
and still compatible with GR at linear order. We will focus on the gravitational potentials 
produced by a "point-like" conserved source. 

The tensor part is described by a massive graviton with mass given by m\. Phenomeno- 
logically, this mass is constraint by cosmological and astrophysical observations (see e.g. 
p9| , |l3|]), and has no impact on the gravitational potentials for point-like sources. Also 
vectors modes do not affect these potentials (for cosmological constraints see pOfl). For our 
purposes only scalar perturbations are relevant. 

Let us briefly review the situation of standard GR in presence of "point like" conserved 
sources, in Minkowski or dS background: 

T 00 = ^, T 0i = Tij = . (7.1) 

In GR, there is no scalar propagating DOF and the gauge invariant potentials are deter- 
mined from the sources as 

$GR = *GR = JfT^ T 00 ■ (7-2) 

Recall that the perturbations are defined with respect to a non-flat metric. Thus, both the 
background and the perturbations play a role in the gravitational dynamics around local 
sources. 

As described in section ^, the generic massive gravity case has two propagating DOF 
in the scalar sector. In this section we are interested in static solutions in the presence of 
static sources. More concretely, we will consider time scales short enough such that we can 
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consider the background metric constant 14 . By inspecting the EOM's in this limit, time 
derivatives can be neglected provided that co, u>H <C k 2 ; loH <C m 2 and luH <C E 2 where 
oj~ 1 is the typical time scale for the variation of the gravitational perturbations and E is 
the energy scale of the sources. 

Once that time derivatives of the two dynamical fields S and r are neglected, and in 
the regime H 2 ,H> \m\ <C mf « A, the EOM can be solved in a straightforward though 
lengthy way. The generalization of the Newtonian potential is the quantity <1> and we get 



712A 2 + mA + no 
d 3 A 3 + d 2 A 2 + diA + d 



(7.3) 



where the nj and di are polynomials in the masses. The physics relevant for the vDVZ 
discontinuity is captured by expanding <I> in powers of 1/A, e.g. A 3> mf . 



00 



u Tnn + v Tu 



M 2 A 



2M 2 A 2 (m 



m% 



+ 0( A3 



u 

V ■ 



[m 4 + m\ (67713 — 4777,4 ~~ 2m|)] , 



a 2 [m 4 , — 2777-2777. 



(7.4) 



Thus, at small distances we get the GR result plus corrections. 15 Also ^, that is important 
for post-Newtonian tests, has the same structure: 



' 00 



Mp A 



2M 2 A 2 



777 



t — 2777n 777 2 



f 00" 



mz 



777o 



7775 



ml 



777$ 



+ 



1 

A3 



(7.5) 



Clearly, no discontinuity is present at small distances provided that m 2 7^ m\ (notice also 



that mi has disappeared from the previous expression). When m\ 



m 



3) 



the previous 



expressions are not valid and a discontinuity is present, as it can be established by noting 
that in the UV the EOM imply 

(7.6) 



2mf* 



which does not hold in GR. 



7.1 Coupling to matter for mi = 

The case m\ = is of particular interest, as in flat space there is no scalar DOF and the 
potential features a correction linear with r, invalidating the linearized approximation at 
large distances. In a curved space the scalar r propagates and the gauge invariant potentials 



14 In this limit the standard Fourier analysis is well suited to analyze the EOM and energy is a conserved 
quantity. It is also clear that if the limit is not singular the results are equivalent to those of Minkowski 
space considered in 

15 The expression 



is valid for distances smaller than the inverse of mass. For distances of the order 



of the inverse of the mass, the appearance of a pole in (7.3) makes the series ill defined. The exact solution 



can be easily found and one can see that the perturbations acquire a Yukawa tail. Thus, this modification 
of Newtonian potential has the desirable feature of keeping the perturbations small at large distances. 
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and can be written as a combination of the source, r and its time derivatives as 
= V GR + a 



2a Hm\m\ r + r' 



$ = i]/ + am 



2AH(m% - ml) 
2 / 2aH{m\ — 3m|) r 



4 / 

mj t 



2 » Afl"(m| - mf) 



(7.7) 



Here we have used the expression for t" obtained from the EOM, namely: 

r" = 2(m '~!!| )g/ (r o - Ml At) + fiK.ff) r + q 2 (m u H) r' , (7.8) 

where 51 2 are functions of the background and the masses, finite in the limit rrij — > 0. 
From these expressions one can study the behavior of potentials in the limit — > 0. 

First, consider the dS background, H' = 0. In this case, the first term in the r.h.s. 
of ( |7.8j ) vanishes. As a result, the only particular solution (vanishing for zero sources) is 
r = 0, and the potentials ( |7.7| ) coincide with those of GR. Remarkably, also the linear term 
appearing in Minkowski ||] is absent in dS. 

A similar situation happens in FRW background: when H' 7^ 0, the first term in the 
r.h.s. dominates in the mi — ► limit, and r remains finite: 

T ~^fF + 0(ra2) - (7 - 9) 

This implies that the corrections to and $ with respect to GR vanish in this limit, and 
there is no vDVZ discontinuity. 

Further insight can be gained by looking at explicit solutions of (|7.8[ ). These can be 
found by assuming a special time dependence of the masses and the scale factor: 

1 mf{r ) ) = a s X i , (7.10) 



where Aj are constants of dimension two. The EOM for r ( |7.8| ) then reduces to 

T + + M|Ay + — + 1)(A2 " A3)A+ 

A 2 [A 4 + i(\2 - 3A 3 ) + A 4 ^(s + 4)]) r = , (7.11) 



,2+s 



where A 2 = \\ + A (A 2 - A3). 

For a general background one can easily find an exact solution of ( |7.11| ) when s = —2. 
The solution that is relevant to us can be written as 16 

Tqq 2 _ 2A 2 A 4 (1 + 21) - (£ + l)Aj + 2l(A| - 3A 2 A 3 ) 

Tp -M2(A-/,2)' " " 2 (^ + i) (A2 _ A3) • 



id 



The general solution is of the form 



+ rp. 



Stability requires |/3(A)| — (2Z + 1) < 0, which at high energies implies i(i + l)A 2 (A2 — A3) oc — ff'A 2 (A2 — A3) > 
0. This condition was readily derived in section 5. J from direct inspection of the Lagrangian. One can also 
check that the solution is stable at any scale if in addition I > —1/2 and A2A 2 {A4 + i(A2 — 3A2 + 2A4)} < 0. 
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In this case, the potentials ( 
/ 2A2A4 

* = %GR + 



arc 



A 2 



2A(A 2 - A 3 ) 



Tp, $ = + A 2 



A2 — 3A3 + A4 
A(A 2 - A 3 ) 



Tp. 



(7.13) 



and one can check that there is no discontinuity in the massless limit. As recalled in 
appendix |b| one can work out the explicit expression for the potentials in position space to 
get (we assume [i 2 > 0) 



* = f 



GR 



<£> = $ 



GR 



1 



1 + 



2A2A4 



A 2 



-fir 



2/x 2 (A 2 -A 3 ) j v 

6A2A3 + 4A2A4 



1) 



2A1 



2/i 2 (A 2 -A 3 ) 




1) 



(7.14) 



where [i can be read from ( 7.12 ). This result differs from the one found in flat space (see 
e.g. in some essential facts: first, instead of the linear correction to <3? that appears 
in Minkowski, we found an exponential function that decays to a constant at large values 
of r (r ^> fi^ 1 ). This is an infrared modification of GR whose magnitude depends on a 
ratio of masses, i.e. it gives finite 0(1) value in the generic m — > limit. Second, also 
is modified at large distances, and the modification decays to a different constant. This 
implies that the modification is not simply a redefinition of Mp. Finally, in FRW t is an 
ordinary propagating DOF (see footnote ^) and, in contrast to the Minkowski case (see 



e.g. [24]), there is no free time-independent function in the solution. At short distance 
both potentials reduce to GR and there is no discontinuity. 

With the above explicit solution one can check that in the dS limit (i — > —1) the 
potentials reduce to GR, because fj, — > 00 and r — » 0, in agreement with the previous 
discussion. On the other hand, also the flat limit £ — * can be safely taken in the last 
expression, but the result is not the Minkowski one. We conclude that the presence of a 
curved background removes the linearly growing term at large distance, or in other words 
regulates the infrared modification of the gravitational force. 

One can also find the exact expression for the potential in the phases where there is 
no scalar DOF, e.g. m v = (see section |6.4.2| ). As we have seen, in this case the kinetic 
term of r is zero and we we can explicitly solve for its EOM for any source obtaining an 
expression similar to ( |7.7| ). Moreover, r will be of the form 



T< 



00 



Ml (A - M 2 ) ' 



M s 



m 



3\2 



H' 



(7.15) 



where q[m,i, a] is an analytic function of the masses, a and their derivatives. The correction 
to Newtonian potential can then be written in the form 



$ = $ 



GR 



V = V 



GR 



1 + 



1 + 



a 777,2 k± 



(ml 



ml)M 2 



a 777n 777^ 



(ml 



mfjM 2 



-Mr 



-Mr 



1 + 



777n 



3m| + 2ml) 



1 



Mre 



-Mr 



1) 



(7.16) 
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where 



,2 2 2^ ml(M 2 Y 

k x = (m| - 3ml + 2m l) + ! 



aHM 2 ' 

Again we see that the presence of a non trivial background gives rise to a modification of 
GR at large distances r ~ M" 1 (pushed to infinity for vanishing masses). 
At short distances, the potential reduces to GR plus corrections: 

* = *GR + ~ 3 ^ + 2 + Q f^_] (717) 

On the other hand taking the dS limit carefully we recover $> = &gr- 
7.2 Coupling to matter for mo = 

For the mo = case, one can express the potentials in a way similar to ( |7.7| ) and ( |7.8[ ), this 



time in terms of E. The resulting expressions turn out to be very complicated, and here we 
will consider explicitly only the "partially massless" case discussed in section |6.2.2| , where 
no DOF is present. In this case, we can write the gravitational potentials as 

.' (3T 00 + Mj,AE)(am| + 2m\H' - 2#(m 2 J) - 2a(3T 00 + 2M 2 AZ)H 2 m 2 



4A 2 H 2 M 2 



3T 00 + MJ.AE 



2A 2 M 2 



* = * GR + a 2 m 2 ( ~Z A ;'J2 ) > (7-18) 



where E is 



M 2 AE = 3T °°( am 4 + 2m l H ' ~ 2 EMl ~ ± aH2m l) (719) 
P a[mi + 2H 2 {2ml-3ml)] + 2[mlH' - H{ml)'}' { ' ' 

The previous two equations indicate that there is no vDVZ discontinuity and we recover 
GR in the massless limit. 

Some care is needed in the special cases when the numerator or the denominator of 



( 7.19 ) vanishes. If the denominator vanishes the theory has an extra gauge invariance (cf. 



section 3.2.2 ). As a result the EMT is coupled consistently only if Too = unless also the 
numerator vanishes. In any case, E can be set to zero by a gauge transformation, and the 
potentials can be read from ( |7.18| ). 

One can readily see from these expressions that the corrections to the Newtonian 
potential simply amount to a linear correction (see appendix [B]) , that invalidates the linear 
approximation at large distance. This modification vanishes for = 0, which also gives 
mi = and H(3ml - m$) = (cf. (ggg)). 



8. Discussion and Conclusions 

In this work we have performed a systematic study of Lorentz breaking massive gravity 
in a FRW background. For the tensor and vector sectors, the analysis is very close to 
that of Minkowski space: both sets of modes acquire independent masses constrained by 
phenomenological bounds. For vector modes, the naive strong coupling scale is similar to 
that of flat space (see also below). 
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The most interesting results are in the scalar sector where generically there are two 
propagating degrees of freedom. For maximally symmetric spaces, the study of the disper- 
sion relations at high energy reveals the appearance of ghost-like instabilities, i.e. instabil- 
ities associated to an infinite volume phase space, that can be cured only by introducing 
a momentum space cutoff. Remarkably, this is not necessarily the case in arbitrary space- 
times: high energy instabilities can be absent in a FRW background with expanding hori- 
zon, i.e. H' < 0, see ( |6. 10| ) and ( |6.12| ). Indeed, at high energies the sign of the determinant 
of the mass matrix A4 (see (|6.8| ) for the definition) is fixed in Minkowski and dS whereas in 
FRW H' enters in the game allowing a region in the parameter space where M is positive 
definite. As a drawback, in the limit H' — > 0, the theory is strongly coupled in the scalar 
sector. 

The scalar sector also features a number of phases with less than two DOF. Gener- 
ically, in a FRW background the phase rriQ = (which includes the FP phase) has one 
scalar DOF. We found the conditions that make the kinetic term of this mode positive def- 
inite, generalizing the Higuchi bound to LB masses in FRW spaces (see ( |6.14| ) and ( |6.15| )). 
Moreover, provided that m\ > m§, high energy instabilities are absent. We also sketched 
the method to avoid instabilities at intermediate momenta (which may be even interesting 
for cosmological perturbations). 

In the presence of curvature there exist situations where the Lagrangian for the scalar 
modes becomes particularly simple as discussed in section 3. 2.1 . In particular the invariance 
under time cliff can be recovered even when mi 7^ 0. More interesting is the case where the 
absence of scalar DOF is due to a residual gauge invariance which is absent in flat space 
(partially massless case). Taking the FP limit in dS, the condition for having residual 
gauge invariance can be solved, and as a result all masses are determined in terms of the 
curvature scale (cf. ( |6.29D ), 

Also interesting is the phase mi = where in general there is again a single propagating 
scalar DOF. For maximally symmetric backgrounds, the EOM for this scalar do not contain 
any gradient term and this mode effectively has no dynamics (zero velocity): it behaves as 
a collective mode. For generic FRW with expanding horizon, the propagating scalar has a 
dispersion relation that can be made free of instabilities. 

We have then analyzed how the Newtonian potentials generated by conserved point-like 
sources are modified. In the general case they agree with GR modulo corrections at scales 
related to the massive gravity scale. A typical form of the gauge invariant gravitational 
potentials is for example (see eqs. ( |7.3|) , ( |7.14j) , (|7.16|) ) 

<D = <D Gfl [l + Cl (e-^ r -l) + C2 /i 2 re-^ r ] , (8.1) 

where the mass scales \i{ are combinations of curvature and mass parameters while q are 
dimensionless combinations. This form is valid also for the 1 f potential. Therefore in the 
massless limit or for scales r <C (aH)^ 1 , m^ 1 the potentials reduce to the GR result, which 
makes these phases potentially very interesting. Comparing for instance to the m\ = 
phase in flat space, where a linearly growing term invalidates perturbation theory at large 



distance [13|, in curved space the potential (|8.1| ) is well behaved at large distance without 
imposing any fine-tuning in the mass parameters. In this sense, the presence of a curved 
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background regularizes many of the peculiarities of Minkowski (also the case m v = 0, 
singular in Minkowski, is regular in FRW). At short distance, the corrections with respect 
to GR in ( |8.lD can be estimated by expanding the exponentials. 

In our analysis we found that some of the propagating states can have small kinetic 
terms (typically proportional to mass or H') giving rise to strongly coupled sectors at very 
low energy. This fact can be relevant for its possible cosmological implications. For instance 
already in the vector sector when m\ ^ 0, perturbation theory is reliable only for H < 
\pm\Mp. This casts serious doubts on the possibility to use massive gravity in standard 
inflation, while keeping small the LB masses. In fact, to trust the linear approximation at 
the standard inflation scale -^inflation — 10 13 GeV one would need m\ > 10 16 eV, and this 
would require a severe fine tuning with respect to the other masses that are constrained by 
various gravitational tests (pulsar, solar system tests) to be much smaller (typically 10~ 21 
eV). On the other hand, the value of ^inflation is very model dependent, the only real upper 
bound comes from BBN, Trh > lOMeV, i.e. H > 10 -16 GeV. This gives the limit mi > 
1CP 30 eV, well below any other gravitational constraints on the masses. The cosmological 
constraints coming from the analysis of the modified gravitational perturbation's dynamics 
are presently under study J2(J] . 

Let us close with a comment on exact solutions. Besides the large distance modifica- 
tions to GR found in this work at linearized level, some modifications have also been found 
in exact (spherically symmetric) solutions of massive gravity [20, 41]. These solutions exist 
in dS space and feature a nonanalytic r 7 term in the gravitational potential. 17 Thus, they 
differ also asymptotically from the linearized gravitational potentials found in this work, 
which may be understood from the presence of long-range instantaneous interactions at 
linearized level. Therefore also for many of the phases analyzed in the present work, one 
may expect important non-linear effects even at large distances. 

Finally, our study suggests that the analysis of perturbations around other nontrivial 
backgrounds may also unveil phases where the perturbations have a stable spectrum, and 
is thus of definite interest. Of main importance would be a dedicated study addressing 
perturbations and their stability in the (exact) gravitational background produced by a 
star. 
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A. The phase m = in the generic case 



In this appendix we present the explicit analysis of the dynamical degrees of freedom of 
the case tuq = discussed in section |6.1[ The EOM for r yield the constraint 



D 



4A (m? - 2ml) H + s \ 2 ( 4A - 3a 2 m 2 ) (am^F 2 - {m\)'H + m\H') 



D = 4A 2 



amf + AH' 



12A 



+ am\ml (2A — 3a 2 m|) 



2am 2 l H 2 - 2(ml)'H + am\m\ + (mf + 2ml) H ' 



(A.l) 



+ 9mf 2am H + a(m|) 2 + 2(ma#' - {mi)' H) 



where m 2 . 



3(mJ 



-m/j 



-m|. Once r is substituted in the action, we get a (quite involved) 



effective action for £ whose kinetic part is written in ( 6.13| ). Requiring the positivity of 
the kinetic energy (no ghost) for large momenta we find the condition ( 6.14| ). As we saw, 
the kinetic term is always positive at small momenta. 

To understand when K, is positive also at intermediate momenta, first notice that K, is 
expressed as a fraction of two polynomials with different roots 18 in the variable A. For the 
fraction to keep its sign those roots must be either at A > or be absent. The numerator is 
a linear polynomial, and setting its root at positive A corresponds to the condition ( 6.15| ). 
The denominator is a second order polynomial that, once the coefficient of A 2 is taken as 
a common factor, has a the zeroth order term 



9a 4 ml 



' a(m\ + 2H 2 mf l ) + 2m\W - 2H{m\) r 
am\ + AH' 



(A.3) 



which, from (|6.14j ) and (|6.15| ) is positive. Thus, the product of the roots of the polynomial 
is positive. Besides, the term proportional to A in the denominator reads 



/ a{m\ + 2H 2 m 2 l ) + 2m\H' - 2H(m\ 



9a 2 (amf + AH') 



12a 2 



am 2 (m 2 — m 2 ) + m 2 H' 
ami + AH' 



(A.4) 



which from ( |6.14D and fl6,15|) is negative definite for m 2 > 0. This finally means that in 
the case m 2 > (required for stability of the vector sector), both roots are positive and 
instabilities in the kinetic term are absent at any scale provided that the inequalities ( 6.1-j ) 
and (|67l5D are satisfied. 

Concerning the potential term, the analysis is more involved. The expansion in A 
provides a useful tool to analyze the absence of instabilities at any momentum scale. The 



3 The two roots coincide only when 



mi(mi — 2ml) 2 [a(m4 + 2H 2 m^) + 2m\H' — //(ml)'] = 0. 



(A.2) 



If the second factor cancels, we find a very simple kinetic term which is always positive. When the last 



factor cancels, the kinetic term is positive at any scale provided that (3.14) holds. Finally, when both terms 
cancel, the constraint (A.l) reduces to r = and the whole action is much simpler. 
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potential term A4 can be written as ( |6,18| ). The absence of gradient instabilities in the 
ultraviolet and infrared regimes, equivalent to the positivity of M in these regimes, requires 



I at large momenta mk > m\ 

< 2 3 

I at small momenta m| > . 

At intermediate scales, the instabilities are Jeans-like. The potential Ly is quartic in A 
which does not allow to find its zero exactly. Nevertheless, imposing that it is free from 
instabilities at high energy scales and at zero momentum we know that it will be positive 
definite at any scale provided that its minima in the interval A 6 (—00, 0] are below zero. 
These minima can be exactly localized as they corresponds to the solutions of 

c + 2dA + 3eA 2 +4(m|-m 2 ; )A 3 = 0. (A.6) 

From the fact that we have at much two minima localized in the interval A £ (—00, 0], and 
yet some extra freedom in the choice of the mass functions, we expect to find a large class 
of lagrangians with a well defined potential (see ||]). 

B. Gravitational Green's Functions 

Once the Newtonian potentials are worked out in momentum space, they can easily be 
Fourier-transformed to the physical position r-space (r = \x\). The potentials <3? and 
found in this work are always of the kind 

= Polynomial^ 11 + • • • ) = ^ gj 

Polynomial(A n + 2 + • • • ) ^ (A - M 2 y 1 ' ] 

where Z{ and Mj are functions of the background and mass parameters. Once the fraction 
has been decomposed in poles, we can use the following correspondence to directly read 
the r dependence: 



1 1 1 



etc, (B.2) 



A r' A 2 A — m? r (A — m 2 ) 2 m 

for suitable choices of integration constants. 
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